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I. INTRODUCTION 

The prediction through the SNela Hubble diagram [1,2] that the universe is at present going 
through on accelerating phase, has been confirmed by wide range of data, from more recent SNela 
data to BAOs and CMBR anisotropies [3-7]. However, the overwhelming abundance of observational 
evidences for cosmic speed up does not match with standard cosmology in the frame work of General 
Relativity (GR). To resolve this paradox, an unexpected new ingredient in the form of a negative 
pressured component has been introduced but it poses difficult questions on its nature and introduces 
further problems hard to be solved. So it is naturally speculated that these observational evidences 
may be the first signal of a breakdown of our understanding of the laws of gravity on cosmological 
scale. As a result several modified gravity theories has been developed of which f(R )-gravity the¬ 
ory [8,9] gets much attention. In this theory the scalar curvature R is replaced by a suitably chosen 
function f(R) in the gravity Lagrangian. 

An equivalent formulation is represented by teleparallelism where instead of curvature, torsion is 
responsible for the gravitational interaction [10-12], This model was first proposed by Einstein for 
unifying electromagnetism and gravity on Weitzenbock non-Riemannian manifold. The Weitzenbock 
connection replaces the Levi-Civita one on the underlying Riemann-Cartan space-time. Thus grav¬ 
itational interaction is not purely geometrical rather, the torsion acts as a force, allowing the inter¬ 
pretation of gravity as a gauge theory of the translation group[13]. Although, conceptually, GR and 
teleparallel gravity are distinct but they yield equivalent dynamics at least at the classical level. 

In analogy with f(R)- gravity, a generalization to teleparallel gravity is obtained by replacing T (the 
torsion scalar ) with a generic function f(T). A particular important consequence is the breakdown of 
the equivalence with the classical GR and hence predicts different dynamics [14-16]. Further, modified 
teleparallel gravity preserves the held equations to be still second order in the held derivatives in 
contrast to fourth order equations in /(R)-gravity. However, this modified theory suffers from the lack 
of Local Lorentz Invariance (LLI), hence all the sixteen components of the vierbien are independent 
and one can not hx six of them by a gauge choice [17]. 

For the past few years, models based on DE interacting with dark matter ( DM ) or any other 
exotic matter components have gained great impetus. Such interacting DE models can successfully 
explain numerous cosmological puzzles namely, phantom crossing, cosmic coincidence and cosmic age 
problem [18-23]. 

The present paper is devoted to the study of dynamics of interacting dark energy in f(T) cosmology. 
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The interaction between DE and DM could be a major issue to be confronted in studying the physics 
of DE. However, due to the nature of these two components remaining unknown, it will not be possible 
to derive the precise form of the interaction from first principles. One has to assume a specific coupling 
from the outset or determine it from phenomenological requirements. Further, in the framework of 
field theory it is natural to consider the inevitable interaction between the dark components. An 
appropriate interaction between DE and DM can provide a mechanism to alleviate the coincidence 
problem. Moreover, complementary observational signatures of the interaction between DE and DM 
have been obtained from the cosmic expansion history by using WMAP, SNIa, BAO and SDSS data 
as well as the growth of cosmic structure. Interestingly, it was disclosed that the late integrated 
Sachs-Wolfe effect has the unique ability to provide insight into the coupling between the dark sectors. 
Further, in view of the continuity equations, the interaction between DE and DM must be a function 
of the energy densities multiplied by a quantity with units of the inverse of time which has the natural 
choice as Hubble parameter. 

Usually one assumes a phenomenological form of interaction between matter and dark energy which 
are the dominant components of the cosmic composition [24-28] today. In fact, the interaction term 
gives the rate of exchange of energy density in the dark sector. 

The motivation for choosing such a complex system is to test of such a weird model could explain 
the overall cosmological evolution. Due to complicated field equations, we study the dynamics in the 
phase space associated to this scenario around both hyperbolic and non-hyperbolic critical points. 
Considering first order perturbations near the critical points, we examine the nature of critical points 
by considering the eigenvalues of first order perturbed matrix. Previously, similar dynamical system 
analysis has been done by us in Dirac-Born-infeld gravity theory [29] and in Brane Scenario [30]. In 
DBI model, critical points are evaluated, analyzed and their stability have been discussed. In brane 
gravity model which is essentially on higher dimensional theory, in addition to the above analysis the 
classical stability of the model itself has been examined. In the present f(T) gravity model, choosing 
the function ’/’ appropriately the evolution equations are converted to an autonomous system by 
suitable choice of the auxiliary variables. In addition to the above course of analysis in the present 
work we shall also examine the stability of the system both classically as well as quantum mechanically. 

We follow the plan : In section II, we present the basic theory of f(T) gravity. In section III, we 
introduce the basic equations and formation of dynamical system. While the critical points of the 
dynamical system and phase space analysis to the corresponding critical points are shown in section 
IV. In section V and VI, we work out the stability analysis of the critical points. The cosmological 
implications at the critical points are discussed in section VII. We conclude and summarize in section 


4 


VIII. Throughout the paper we use natural units ( 8nG = - 4 ?— = h = c = 1 ). 

m PL 

II. /(T)- GRAVITY THEORY 

A possible way to modify gravity beyond general relativity is to consider the Weitzenbock con¬ 
nection, which has no curvature but torsion, rather than the curvature defined by the Levi-Civita 
connection. Such an approach is termed as ’teleparallelism’[12,31-33], which was also considered by 
Einstein himself[10,ll]. However, in the context of the present accelerated expansion of the universe, 
the teleparallel Lagrangian density described by the torsion scalar T has been extended to a function 
of T ( f(T) ) [34,35], equivalent to the concept of f(R) gravity. 

In this new modified gravity theory ( called f(T)- gravity ), the gravity is no longer caused by curved 
space-time but torsion and moreover, the field equations are only second order unlike the fourth order 
equations in the f(R) theory. 

In f(T) gravity theory the action is written as 

I= W J ^[M(T + /(T)) + J L m ] (1) 

Where T is the torsion scalar, f(T) is a differentiable function of the torsion, L m corresponds to 
the matter Lagrangian, |e| = det(e^) = sj—g and k 2 = 8nG. The torsion scalar is defined as 

T = S p ^TP ^ (2) 

with 

S p ^ = X -(K^ p + 5^T 9u e - SpT 9 ^ g ), (3) 

P = -^(T^ p-T^ p-Tp^) (4) 

T x , w =T u, - r pu = - d v e*) (5) 

In teleparallelism, orthogonal tetrad components eA(x^) are used as dynamical objects and are 
cosidered as an orthonormal basis for the tangent space at each point of the manifold. So we 
have 


&A&B = VAB = diag( 1, —1, —1, —1) 
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. Here each vector e A can be described by its components , /r =0,1,2,3 (A =0,1,2,3) in a coordinate 
basis i.e. e A = ^ A d^- Here capital letters refer to the tangent space while greek indices label coordinates 
on the manifold. The metric tensor is obtained from the dual vierbein as 

9iny{x) = rj AB e^{x)e^(x) (6) 

In the above T x )W is the curvature less Weitzenbock connection and it encompasses all the information 
about the gravitational field. The contorsion tensor K^ u p gives the difference between Weitzenbock 
and Levi-Civita connections. 

It is straight forward to show that this equation of motion reduces to Einstein gravity when 
f(T) = 0. Indeed, it is the equivalency between the teleparallel theory and Einstein gravity [12]. The 
theory has been found to address the issue of cosmic acceleration in the early and late evolution of 
universe [36] but this crucially depends on the choice of f(T). For instance, exponential or power-law 
form of f(T) can not lead to phantom crossing [37]. But subsequently, it has been shown that crossing 
of the phantom divide line for the effective equation of state of two specific models in f(T) gravity 
is possible [38] and the best fit results suggest that the observations favour a crossing of the phantom 
barrier. Then reconstruction of f(T) models has been reported in [39,40] while a detailed cosmological 
analysis is performed in [41,42] and thermodynamics of f(T) cosmology [43] including the generalized 
2nd law of thermodynamics has been investigated. 

Further, it should be noted that a constant f(T) acts like as a Cosmological Constant while / 
linear in T (i.e fr = constant ) is simply a redefinition of Newton’s Constant G. A desirable choice 
of f(T) is such that at high redshift general relativity should hold so as to agree with primordial 
nucleosynthesis and cosmic microwave background constraints i.e f /T —>• 0 at early times (a <C 1 ) 
[44]. On the other hand, in the asymptotic future there will be de Sitter state. A simplest choice of 
/ is the polynomial form [44] 


f(T) = P(-T) n (7) 

where j3 is arbitrary constant. 

In particular for n = ^ the model gives the same expansion history as DGP gravity [44,45] and 
hence f(T) gravity can be considered to be related to higher dimension theories. Note that the 
restriction 'n <C V gives /(T) a viable model compared to current observed data set and also puts 
the rescale factor fc to Newton’s constant to be small. The effective DE equation of state varies from 
u) = — 1 + n in the past to uj = — 1 in the future. 
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Throughout the work we consider the flat homogeneous and isotropic FRW universe with the metric 


ds 2 = dt 2 — a{t) 2 'y^ jd'j 


,i \2 


( 8 ) 


where t is cosmic time. For this metric 


= diag(l,a(t),a(t),a(t)), 

where a(t ) is cosmological scale factor. By combining with (5), (3) and (4) one obtains 


T = S ppv T pllv = -6 H 2 , 


where H = - is the Hubble parameter and dot stands for derivative with respect to t. 


III. BASIC EQUATIONS IN f(T) COSMOLOGY AND FORMATION OF DYNAMICAL 

SYSTEM : 

In f(T) gravity theory, we consider flat, homogeneous and isotropic FRW space-time as the model 
of our universe and the matter is chosen as dark matter in the form of dust (having energy density 
p m ) interacting with the dark energy which is chosen as a minimally coupled scalar field cj) having self 
interacting potential V(4>). So the energy density and pressure corresponding to this scalar field are 

P<t> = + V(0),P<i> = 7^ 2 - V(4>) (9) 

Thus the modified Friedmann equations in /(T)- gravity are[14,15] 


H 2 


(2f T + l) [ 3^ + Pm) 6 ] 


( 10 ) 


^ _ 1 r Pm + P4> + P<t> 1 

2 1 + /t + 2T fxT 

The energy balance equations for the individual dark components are 


Pm T 3 H p m — () 


( 11 ) 


( 12 ) 


and 
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p<t> + 377(1 + uj^p^ — —Q (13) 

Where the interaction term Q corresponds to energy exchange between dark energy and dark matter. 
The positivity of Q indicates a transfer of energy from dark energy to dark matter. This is required 
to alleviate the coincidence problem and is compatible with the 2nd. law of thermodynamics. For 
the time being Q is unspecified, only it is assumed that Q does not change sign during the cosmic 
evolution. Due to the unknown nature of the two dark components (DM and DE), the precise form 
of the interaction can not be determined from the outset or from phenomenological requirements. 
However, it is speculated that the interaction term may lead to a major issue to be confronted in 
studying the physics of DE. Also, from the view point of field theory, interaction between the dark 
components appears naturally. Moreover, an appropriate interaction between the DE and DM can 
provide a mechanism to alleviate the coincidence problem. Further, in view of the continuity equations 
( i.e. equations (12) and (13) ) the interaction between DE and DM should be a function of energy 
densities multiplied by a factor having dimension inverse of time and Hubble parameter is a natural 
choice for it. Thus phenomenologically, Q can be chosen as (i) Q = Q(Hp m ) , (ii) Q = Q^Hp^), (iii) 
Q = Q[H(p,p + p m )\ or more generally (iv) Q = Q(Hp^, Hp m ). In the present work, for simplicity Q 
is chosen as Q = aHp m , where the coupling parameter 'a' is assumed to be small. 

Now, using (9) in the continuity equation (13) the evolution of the scalar field is given by 


<j) + 3 Hcj) + 


dV{</>) 

dc(> 


(the modified Klein-Gordon equations) 


Q 


(14) 


Due to complicated form of the evolution equations, it is not possible to have analytic solution so 
to have a qualitative idea about the cosmological behavior we shall put the evolution equations into 
an autonomous dynamical system. For this, we introduce the new variables 


x = 


_ y/m n _ Pm 
> y — /— TT > ^ L m — 


y/6 H' y V3H 


3 H 2 


(15) 


which are normalized over Hubble scale. is the density parameter for DM. As a result, the 
evolution equations reduce to the following autonomous system of ordinary differential equations ( 
after some algebra ) 







dx 

dN 
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= (1 + x 2 - y 2 )(^x - ^-) + (a - 3)x - y^y 2 A 

—— = y[\f^-x\ + -(1 + x 2 — y 2 )\ 
dN yL V 2 2 v ;J 

= n m [a + 3(x 2 -y 2 )} (16) 

where we choose \/^j = A, a constant such that V(<t>) = e x ^ and the independent variable is chosen 
as N = In a, which is called the e-folding parameter. The system of equations in (16) is analyzed by 
first equating them to zero to obtain the critical points. Next we perturb equations up to first order 
about the critical points and check their stability. 


In deriving the above autonomous system of ODE we choose f(T) = j3\J—T (i.e n = ^ in equation 
(7) ) and the interaction is chosen in the form Q = aHp m . As mentioned earlier, the positive coupling 
is very reassuring in view of coincidence problem. However, models with negative coupling parameter 
(indicates decay of DM into DE ) allow for the possibility that there is no DE field in the very early 
universe and that DE ’Condenses’ as result of the slow decay of DM [46]. Also, it has been shown 
[47] that the coupling parameter is weakly constrained to negative values by Planck measurements. 
Further, the negative coupling can not be counted to resolve the tension between the Planck and 
HST measurements of the Hubble parameter [47]. Although, the negative coupling does not help to 
alleviate the coincidence problem, it appears in the observed data fittings that models with negative 
coupling shows most significant departure from zero coupling. 

Further, using the normalized variables in the first modified Friedmann equation, we obtain the 
density parameter for the dark matter as 


— 1 x y 


(17) 


Due to the energy condition 0 < Q m < 1 , so for fix , (x,y) lies on the circle x 2 + y 2 = 1 — . In 

fact, the phase space (x, y, D m ) of the autonomous system (16) forms a paraboloid ( x 2 + y 2 + VL m = 1 
) bounded by Q m = 0 and = 1. Hence the phase space is finite . 

The cosmological parameters related to the scalar field namely the equation of state parameter 
and the density parameter can be expressed by newly defined variables as 

_p < f,_x 2 -y 2 _ Ptf, _ 2 2 

bJ(f> — | _ 9 ; ^ ~0 0 rJ -9 % T y 


P <t> X z + y z 


3 H 2 


( 18 ) 
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and 


P 

UTot = - = 


p<t> 


x 2 — y 2 


P P<t> + Pm x 2 + y 2 + tt r 


2 2 
= x -y 


since x 2 + y 2 + = 1- Also, the deceleration parameter has the explicit form : 

q = = — 1 + 2 ^ x2 ^ m ) 


and the equation of state parameter of the equivalent geometric matter is 


(19) 


( 20 ) 


U! 
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(g~ 2 ) 
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IV. CRITICAL POINTS AND PHASE -SPACE ANALYSIS : 


( 21 ) 


We analyze the stability of the corresponding dynamical system about the critical points. We shall 
plot the phase and evolutionary diagrams accordingly. For this reason we must find the critical points 
of the system of equations (16) forming an autonomous dynamical system and then we linearize the 
system near the critical points. The system of equations (16) has the following six critical points for 
the positive coupling parameter a i.e. decaying of DE into matter : 


• I. Critical points: P\ . p> = (±1,0,0), 

• II. Critical Points : ± 3 , 1-4 = (— ±\J l — 0) 

. III. Critical points : P 5 ,P 6 = (Jjg, ±/f + ^(1 - ^)). 

where A = jr = constant , V' = 

On the other hand, for negative coupling parameter we have the following two critical points 

• IV. Critical points : ± 7 , Pg = (± \J~ f, 0,1 + (|) where a € (—3, 0) . 

These critical points and the relevant physical parameters at those points are shown in table I and 
table II. 

From table I, we see that the critical points Pi and P 2 always exist ( for all a, A ) while P 3 and P 4 
exist only for A 2 < 6 . Also, the four critical points Pi, P 2 , P 3 and P 4 represent only the DE compo¬ 
nents ( DM is absent ). The critical points P 5 and Pg exist for 2aA 2 + (a — 3 ) 2 > 0 and correspond 
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TABLE I. Tabic shows the location of the critical points and the values of the relevant physical parameters at 
those points for positive a . 


Pi 

X 

y 



^Tot 


q 

U}g 

Pi 

1 

0 

0 

1 

1 

1 

2 

0 

p 2 

-1 

0 

0 

1 

1 

1 

2 

0 

p 3 

A 

a/6 

Aw 

0 

T—1 

A- 1 

3 

1 

-1 +Ai 

x ^ 2 

-1+ V 
T 6 

p 4 

A 

V6 

7 V 1 - 1 

0 

T—1 

Ai_i 
3 1 

1 

\2 

1 ^ 2 

\2 

— 1 4- — 

p r 

a—3 

, lot , (a-3) 2 

3—a /1 3—ol \ 

3 V 1 A 2 ) 

aA 2 

a 

a. 

(a-3) 2 

1—a 

(3+a) 

-*5 

vdiA 

V 3 1 6A 2 

(a-3) 2 +aA 2 

3 

3 

1 3A 2 

2 

6 

p 

a—3 

Jet | 0-3) 2 

3—a.(-\ 3 —ex \ 

3 ->?- > 

aA 2 

a 

a. 

(a-3) 2 

1—a 

(3+a) 

r 6 

x/6A 

V 3 H 6A^ 

(a—3) 2 +aA 2 

3 

3 

H 3P 

2 
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TABLE II. Table shows the location of the critical points and the values of the relevant physical parameters 
at those points when — 3 < a < 0 . 


P 

X 

y 



WTot 

£l(j) 

q 

Ulg 

P7 

a / — 

0 

x ^ 3 

1 

ol 

3 

a 

3 

|(1-3a) 

(1+a) 

2 

p 8 

^1 - a 

0 

1 + 2 
x ^ 3 

1 

a. 

3 

a 

3 

1(1-3a) 

(l+a) 

2 


a combination of DM and DE with the ratio of two energy densities r = . For 

all critical points DE behave like a perfect fluid. Further, in table II we see that both the critical 
points P-j and P$ are combination of DM and DE with the ratio of density parameter = —1 — A 
The equivalent geometric matter behaves as dust for the critical points Pi and P 2 , it behaves as 
quintessence matter for the critical points P 3 and P 4 while the geometric matter acts as phantom fluid 
for the critical points P 5 and Pq. 


V. STABILITY ANALYSIS 

We shall now discuss the stability of the critical points ( presented in Table I and the table II ) of 
the autonomous system (16), considering first order perturbations near the critical points. To examine 
the nature of critical points, one has to study the eigenvalues of the first order perturbation matrix 
which have been presented in Table III and table IV below : 


From table III, for A = —\/ 6 , the critical point Pi is non-hyperbolic in nature otherwise hyperbolic. 
On the other hand, for A = JE, the point P 2 is non-hyperbolic in nature otherwise hyperbolic. 
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Vector field and Trajectories of the autonomous system (16) 



-2 -1.5 -1 


FIG. 1. Projections of the phase trajectories onto the plane ( x, y, = 0 ) of 
interacting DE in f(T) model for the choices of (a = 0.001 , A = 2.7) 


Vector field and Trajectories of the autonomous system (16) 



FIG. 2. Phase portrait of the system (16) for the choices of a = 0.001, A = —2.7. 
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FIG. 3. Figure shows the phase portrait of the 
system (16) for the choices of a = 0.001, A = 1.1 
in the phase plane (x,y,f 2 m = 0 ). 


FIG. 4. Phase portrait of the system (16) for the choices 
of a = 0.001 , A = -1.1 


Vector field and Trajectories of the autonomous system (16) 



FIG. 5. Phase portrait of system of equations (16) 
for the choices of a = 0.001, A = 1.9 in the phase 
plane ( x, y, ft m = 0 ) . 


Vector field and Trajectories of the autonomous system (16) 



FIG. 6 . Phase portrait for a = 0.001, and A = 2.1 


The critical point P\ is saddle node for A < — y/E ( see figure 2 . ) and P2 is saddle node for A > y/E ( 
see figure 1 . ) while P\ is unstable node for A > —y/E ( fig: 1 . ) and P2 is unstable node for A < y/E 
( fig: 2. ) There are no acceleration phase of the universe near the critical points Pi and P 2 . Also, 
these two critical points represent solutions without matter part of the universe. They are only the 
kinetic energy dominated solutions. 

On the other hand, the critical points P 3 and P 4 are same in all respect and are hyperbolic in 
nature. These two critical points are stable solutions if 0 < a < 3 — A 2 . Also, they behave like 


































Vector field and Trajectories of the autonomous system (16) 



FIG. 7. Projection of phase trajectories onto a given plane ( x,y ,fi m = 0 ) here, we 
choose fl m = 0 . Note that while the point Pi : (x,y)=(l,0) and the point P 2 : (x,y) 
= (-1,0) are unstable node i.e. source for (a = 0.001, A = 0.5), the points P 3 and 
P 4 (see in table I) are the stable solutions i.e. attractor solutions in the given phase 
space. Observe that all the trajectories in the phase space always emerge from the 
points Pl and P 2 while all the trajectories enter into the points P 3 and P 4 


Vector field and Trajectories of the autonomous system (16) 



Vector field and Trajectories of the autonomous system (16) 



FIG. 8 . Phase portrait of system of equations (16) 
for the choices of a = 0.001, A = 0.5 with grid 
difference 0.52 


FIG. 9. Phase portrait for a = 0.0012, and A = 0. 
with grid difference 0.5 
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TABLE III. 


Eigenvalues of the linearized matrix for the critical points of the autonomous system (16) when a 


is positive and A = [(—| + f + \ ' * A2 3 ^ 


^) 2 + (4a 2 A 4 - (a - 3) 4 )(A 2 + a - 3) ( 5 ^ - 3)]* 


P 

Ai 

^2 

CO 

Pi 

3 H- OL 

3 + \/i A 

3 + a 

P2 

3 H - Oi 

3 — \/f A 

3 + a 

P 3 

a — 3 + A 2 

A_ ‘X 

2 0 

a — 3 + A 2 

Pa 

a — 3 + A 2 

Ai -3 

2 6 

a — 3 + A 2 

P 5 

0 

1 [ A 3 3c>! aX 2 1 3(a— 3) 2 i 

2 c 1 2 2 a—3 ' 4A 2 J 

1 r A 3 3a aA 2 , 3(a—3 ) 2 1 

2 1 2 2 a—3 ' 4A 2 J 

Pq 

0 

1 r A 3 3a aA 2 3(a-3) 2 ! 

2 v 1 2 2 a- 3 4 4A 2 J 

It A 3 3a aA 2 , 3(a-3 ) 2 1 

2 L 21 2 2 a—3 4 4A 2 J 


TABLE IV. Table shows the eigenvalues of the linearized matrix for the critical points of (16) when a G (—3,0) 


P 

Ai 

^2 

A 3 

P? 

—a — 3 

f +1 + 

0 

P 8 

—a — 3 

a , 3 \ / a 

2 ' 2 A V 2 

0 


attractor solutions for A 2 < 3 ( i.e. — \/3 < A < y/3 ). But when 3 < A 2 < 6 the two critical points P 3 
and P 4 behave like saddle node ( these are shown in the phase portrait ). There exists an accelerating 
phase of the universe near these critical points P 3 and P 4 for A 2 < 2 . 

Pq and Pq are also same in all respect. They are combination of both DM and DE. These solutions 
are dominated by DE components only for a+A 2 = 3. Pq and Pq correspond to solutions in accelerating 
or decelerating phase of universe according as a >< 1 . 

Now, in the table IV, for negative coupling parameter, we get two critical points namely TV and 
Tg which are non-hyperbolic in nature. The solutions are saddle like ( because one eigenvalue is 
always positive in both of the cases ) due to instability in the eigen direction associated with positive 
eigenvalue and the stability of an eigen direction associated to a negative eigenvalue. These solutions 
are dominated by both DM and DE components. There is no accelerating phase near these critical 
points. 


TABLE V. Stability criteria of the model where B = 6a; 2 + 2\^6xy 2 X + afl m > 0. 


c s 2 

For classical stability (C 2 > 0) 

For quantum stability 

(Px > 0 , p x + 2xp xx > 0 ) 

1 1 VExy 2 X 

6a; 2 + 2^6xy 2 X + afl m > 0 

6a; 2 + 2\/6xy 2 X + afl m > 0 and 

1 + ^f[(6x + V6Xy 2 )^ + 2^6A xy% + f ] 

< 9xH(2x 2 + f2 m ) 

' r 3a 2 + f n m. 
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Thus, our main results can be summarized as follows ( see tables I, III, II and IV ): 

• For all A and a > 0, the kinetic energy dominated solutions ( points P\ and P 2 in table I ) are 
always unstable points in phase space. In particular, P\ and P 2 represent source or repealer ( unstable 
nodes ) in phase space for A E (— V 6 , \/ 6 ) and are always decelerating (q = 2) (see figures 3, 4 and 7 ). 

• For 0<a<3 — A 2 , A 2 <3( i.e. —a/ 3 < A < \/3 ); the scalar field dominated solutions ( points 
P 3 and P 4 in table I ) are the late time attractor ( see figures 3 , 4 and 7 ), while for 3 < A 2 < 6 the 
solutions are saddle node ( see figures 5 and 6 ). Also, for A 2 < 2 , there exists an accelerating phase 
of the universe near the critical points P 3 and P 4 . 

• For all A and a > 0 two critical points P 5 and Pq ( see table I ) are non-hyperbolic and are com¬ 
bination of DE and DM components. These are the solutions accelerating or decelerating completely 
depend on the coupling parameter a. By the linear stability theory, it is not easy to get the stability 
criteria of these critical points. 

• For a € (—3,0), the two non-hyperbolic critical points P 7 and and Ps ( see in table II ) can 
be obtained which are saddle points and are both combination of DE and DM. They are always 
decelerating ( q = ^(1 — 3a), see table II ). 

• In the above figures 1-7 the variables x and y are chosen in the range [-2,2] with grid points at a 
difference of 0.5. The model parameters a and A are chosen as a = 0.001 and A = 0.5. To examine, 
whether the present model depends continuously on data we have drawn figures 8 and 9. In fig. 8 
we have changed the grid difference to 0.52 while in fig. 9 we have changed the model parameters as 
a = 0.0012 and A = 0.52. It is found that both the figures 8 and 9 do not differ significantly from 
figure 7. So we may conclude that the model depends continuously on data. 

VI. EQUILIBRIUM POINTS AND STABILITY CRITERIA 

In the present three dimensional autonomous system the local stability criteria of an equilibrium 
point is characterized by the eigenvalues of the perturbation matrix ( presented in table III and table 
IV for different values of coupling parameter a ) and then discussion about local stability is presented 
. We shall now investigate the classical as well as quantum stability of the model. 

In cosmological perturbation, sound speed ( C s ) has a crucial role in characterizing classical 
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TABLE VI. Condition for stability at each equilibrium point. 


Pl 

X 

y 


local stability 

classical stability 

quantum stability 

Pi 

1 

0 

0 

unstable 

stable 

stable if H > jh 

P2 

-1 

0 

0 

unstable 

stable 

stable if If < — -h 

P 3 

A 

\/6 

A-f 

0 

stable (sink) if 

0 < a < 3 — A 2 and A 2 < 3 

stable if A 2 > 3 

stable if A 2 > 3 and 
X 3 H <~\Jl 

Pa 

A 

V6 

-A-t 

0 

stable (sink) if 

0 < a < 3 — A 2 and A 2 < 3 

stable if A 2 > 3 

stable if A 2 > 3 and 

A 

P 5 

a:-3 
V6A 

/a , (a-3) 2 

V 3 + 6 A 2 

^(1-3^) 

linear stability fails 

stable if a > 3 

stable if a > 3 and 

Pe 

a-3 

V6A 

la , (a-3) 2 

V 3 + 6 A 2 

3—af-i 3—a \ 

3 T 5 "" > 

linear stability fails 

stable if a > 3 

stable if a > 3 and 

(«- 3) 2 # < 

Pi 

_ a 

0 

1 + - 

unstable (saddle) 

unstable 

unstable 

Ps 

/ oT 

V 3 

0 

1 + - 
1 ^ 3 

unstable (saddle) 

unstable 

unstable 


stability. In fact, C 2 appears as a coefficient of the term ^ ( k is the comoving momentum and 
’a’ is the usual scale factor ) and classical fluctuations may be considered to be stable when C 2 is 
positive. On the other hand, for quantum instabilities at UV scale we decompose the scalar field into 
a homogeneous part ( 4>o ) and a fluctuation as 


cj)(x,t) = 4> 0 (t) + 54>(x,t) (22) 

Then by expanding the pressure p(x, 4>), up to second order in 5(j), the Hamiltonian for the fluctuations 
takes the form [48,49] 


H — ( p x T ‘^xpxx') 


m 2 

2 


W ) 2 

Px 


p<j><t> 


m 2 

2 


where suffix stands for differentiation with respect to the corresponding variable [50]. 
For positive definiteness of the Hamiltonian we must have 


(23) 


Px H - ^^Pxx ^ 0 iPx ^ 0, P(j)(j) ^ 0 


(24) 
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where the first two inequalities are related to quantum stability. 

Further, C 2 should be less than unity. Because otherwise it is possible to send signals along space¬ 
like world lines and this would open a Pandora’s box of classical time travel Paradoxes and this also 
violates the unitarity of quantum theory. Further, if C 2 > 1 then from the effective field theory 
the fluctuations 5(j) propagate faster than the speed of light and therefore there will be violation of 
Causality. In the present cosmological scenario we have 


ci = \ + 


\/Qxy 2 X 
3x 2 + j Q r 


so for classical stability ( assuming a > 0 


(25) 


6 x 2 + 2\[§xy 2 X + aO m > 0 (26) 

while for quantum stability the restrictions are 


6 x 2 + 2 \/6xy 2 X + a£l m > 0 


and 


1 + -^-[(62 + + 2 \/ 6 Ax 2 /\^ $xH(2x 2 + tt m ) 


where B = 6 x 2 + 2\/6xy 2 X + aQ m which is positive, H is the Hubble parameter and the values of 
^ and are given in the system of equations (16). 

We have shown both classical and quantum stability criteria of the model in table V. We shall now 
discuss about the criteria for the model stability at the equilibrium points ( presented in table I and 
table II ) when x,y and take the corresponding values at the equilibrium points. From the tables 
I and III we see that the equilibrium points Pi and P 2 are not locally stable( saddle node or unstable 
node ), they behave like source for A 6 (— \/ 6 , V&) and from the above model stability analysis they 
correspond to classical stability only but conditional quantum stable i.e Pi corresponds to quantum 
stability if H > yg while for H < — yg, the critical point P 2 corresponds to quantum stability. As the 
equilibrium points P 3 and P 4 are same in all respect so they correspond to classical stability of the 
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model if A 2 > 3 while for quantum stability the restrictions are A 2 > 3 and A 3 H < — y' | ( see table 
VI ). From local stability analysis, we see that the critical points Pq and P4 ( table I ) are stable if 
A 2 < 3. From the above model stability analysis we see that the two equilibrium points P5 and Pq are 
classical stable if a > 3 where as they correspond to quantum stability if a > 3 and (a — 3) 2 y < 

( presented in table VI ). Finally, from the tables II and IV, we see that two equilibrium points P7 
and Ps are not locally stable (saddle) and from the above analysis we see that they are not classical 
as well as not quantum stable. The corresponding condition for stability at each equilibrium point 
are presented in table VI. 


VII. EQUILIBRIUM POINTS AND COSMOLOGICAL IMPLICATIONS 

From the above phase space analysis of the f(T ) model we shall now discuss about the cosmological 
behavior of the model at the equilibrium points. From table I, we find that Pi and P2 indicate the 
universe is completely dominated by kinetic energy of the scalar field and late time acceleration is not 
possible. They correspond to flat non-accelerating unstable universe without matter part and so they 
are not of much interest in the present context. 

On the other hand, the points P3 and P4 are interesting from the cosmological point of view (for 
A 2 < 2 ). There exists an accelerated expansion of the universe near P3 and P4 which correspond to 
DE model in the quintessence era and the critical points are stable as well as the model is locally 
stable. It should be mentioned that these critical points correspond to cosmological solutions which 
describe the recently observed late time acceleration of the universe ( details are shown in figures 
). Thus, our model predicts the correct evolutionary scheme for DE density in the regime of this 
quintessence scalar field. However, for the restriction 2 < A 2 < 3 the points P3 and P4 are stable from 
local analysis and at these points the potential energy dominates over the kinetic part and the scalar 
field behaves as exotic fluid. The equilibrium points P5 and Pq are the potential energy dominated 
non hyperbolic solutions so we cannot investigate the local stability criteria of the system by the 
linear stability theory. They are the combination of both DM and DE. Accelerating or decelerating 
phase of the universe near the points P5 and Pq are fully depended on the coupling parameter of the 
interaction term of DM and DE. The equilibrium points P-j and Pg are same in all respect for the 
negative coupling parameter. These are the kinetic energy dominated solutions of the scalar field and 
there exists an non-accelerating universe near P7 and P$. 
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VIII. DISCUSSION AND CONCLUDING REMARKS 

In this paper we have performed a dynamical system analysis of a complicated model based on an 
alternative theory of gravity- the so called f(T) gravity. The matter is chosen as dark species - dark 
matter and dark energy which are represented by dust and a scalar field with an exponential potential 
respectively. We have investigated the stability and phase space description of interacting dark energy 
in f(T ) cosmology by introducing a simple interaction between dark energy and the matter content of 
the universe. We have written the general dynamical system equations from basic evolution equations 
by the suitable variables which are normalized over Hubble scale and have found two sets of critical 
points for positive and negative coupling parameter ( see table I and II ). In the context of f(T) grav¬ 
ity the negative coupling parameter shows that the energy is being transferred from matter to dark 
energy. We get a finite phase space of the system of evolution equations ( autonomous system ) and 
phase space forms a paraboloid bounded by Q m = 0 and f l m = 1. From the phase space analysis and 
stability analysis of the critical points ( presented in tables I and II ) we find the stable and unstable 
solutions ( see figures 1-7 ) for the different values of parameters involved. Also the accelerating or 
decelerating behaviour of the universe near both the hyperbolic and non-hyperbolic critical points 
have been given in the last section. It is worth noting that the recently observed accelerated expansion 
of the universe near the equilibrium points P 3 and P 4 can be seen from the autonomous system (16) 
in the quitessence scalar field. These support the recent evolutionary scheme of the universe. 

Further, we have investigated the classical as well as quantum stability of the model. Note that 
these two types of stability are not interrelated because the stability of the critical point is related 
to the perturbations 5x , 5y and <5fi m ( the corresponding variables of critical point ). On the other 
hand, the classical stability of the model is connected to the perturbations 5p (and depends on the 
conditions Cg > 0 ) while the quantum stability is related to the perturbations 5(j) and the conditions 
take the form of inequalities (24). Thus the critical points can be classified into three categories namely 

(i) unstable points at which the model is stable . 

(ii) stable points at which model is unstable and 

(iii) stable points with stable ( both classical and quantum ) model. 

From the table I, the equilibrium points P 3 and P 4 are interesting in the present context, they 
purely describe the scalar field dominated solutions in the f(T) model and describe the late time 
acceleration. So, we conclude that although the equilibrium points P 3 , P 4 correspond to stable model 
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configuration but locally they are not stable equilibrium points. However, it should be noted that the 
dimensionless variables introduced in equation (15) are not suitable for a complete analysis. Though, 
the analysis of the static solutions that correspond to equilibrium points at infinity could be analyzed 
by Poincare projection technique, but for bouncing scenario the expansion rate may pass through zero 
making the state space non-compact. Therefore, for future work, it is desirable to choose appropriate 
compact variables to care of expansions, collapses, static solutions and bounces as well. 
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